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Abstract 

We study a stochastic model for the coevolution of a process of opinion formation in 
a population of agents and the network which underlies their interaction. Interaction 
links can break when agents fail to reach an opinion agreement. The structure of the 
network and the distribution of opinions over the population evolve towards a state 
where the population is divided into disconnected communities whose agents share 
the same opinion. The statistical properties of this final state vary considerably as 
the model parameters are changed. Community sizes and their internal connectivity 
are the quantities used to characterize such variations. 
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In a complex system, the spontaneous emergence of collective non-equilibrium 
behaviour, such as coherent spatiotemporal structures or synchronized dynam- 
ics, is driven by mechanisms which involve both the interactions between the 
system's components and external influxes [1,2]. In segregation phenomena, a 
form of self-organization well known to physicists and chemists, an ensemble 
of interacting elements becomes split into subensembles whose components 
share certain individual states. Segregation occurs also in biological and so- 
cial systems [3,4,5], where it plays a crucial role in sustaining diversity at 
many levels -cellular, functional, organizational, ecological, cultural. Though 
it is usually associated with phase separation in space, segregation not always 
takes place in the spatial domain. In human societies, for instance, two or 
more subpopulations or communities may coexist in the same geographical 
region and, yet, exhibit mutually excluding cultural traits [6]. With respect to 
those traits -which may involve religious beliefs, professional or generational 
jargons, artistic inclinations, etc.- such communities are scarcely interacting, 
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and can be effectively considered as segregated from each other. A key mecha- 
nism in this kind of social segregation is the feedback between the construction 
of agreement within a community and the enhancement of distinctions with 
other communities: specific traits become better established as differences be- 
tween communities develop and grow. 

In the last years, physicists have been increasingly interested in the dynamical 
and statistical modelling of complex systems of biological and sociological 
inspiration as populations of agents whose interaction patterns are described 
through graphs, or networks. Much attention has been paid to dynamical 
processes defined over quenched networks [2,7], and to network growth -both 
purely stochastic [8] or driven by adaptive mechanisms [9]- with emphasis in 
the statistical properties of the resulting patterns. On the other hand, the 
possible coevolution of the network structure and the dynamics taking place 
over them seems to have been less studied (see, however, Ref. [10]). This 
kind of coevolution is at the basis of the feedback phenomenon pointed out 
above, where the formation of internally homogeneous subpopulations and the 
weakening of their mutual interactions enhance each other. 

The aim of this letter is to present a very simple model for coevolution of 
a population of agents and their interaction network. The agents' dynamics 
is based on an elementary model of opinion spreading [11]. The population 
starts in a situation where every agent is able to interact with any other 
agent, and evolves towards a segregated state with disconnected communities. 
Interactions between agents with similar or different opinions are respectively 
favoured or penalized. In spite of the simplicity of the evolution rules, the 
population can reach a variety of social patterns, which range from splitting 
into two communities of similar sizes and opposite opinions, to a single large 
community with homogeneous opinion. Typical connectivities and opinion dis- 
tributions, resulting from the combined evolution of the network and the agent 
states, change considerably with the control parameters. 

The system consists of N agents situated at the nodes of a network. Initially, 
the network is fully connected, so that any pair of agents can potentially 
interact. The individual state of each agent is assigned at random one of 
two possible values, say +1 or —1, with equal probability. Individual states 
represent the agents' opinions, which may eventually change by interaction 
between agents connected by the network. 

The evolution runs as follows. At each step, a pair of connected agents is 
chosen at random from the whole population. If both agents have the same 
opinion, nothing happens. Otherwise, with probability pi either agent adopts 
the other agent's opinion, so that the two opinions become identical. With 
the complementary probability 1 — pi, opinions are not changed. In this case, 
however, the link between both agents is broken with probability p 2 , and the 
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interaction network looses one of its edges. These rules are successively applied 
until no further changes are possible. Since network edges are irreversibly 
deleted by the evolution, a frozen final state is eventually reached where, 
generally, the network is split into disconnected subsets. Within each of these 
communities all agents share the same opinion. 

The frequencies of the individual events that drive the dynamics depend on 
the probabilities pi and p 2 . With respect to the evolution of the system, these 
probabilities are independent control parameters. The statistical properties of 
the final state, however, are completely determined by the relative frequencies 
of the different processes that effectively change the state of the system. In 
other words, they depend on pi and p 2 through certain algebraic combinations 
only. To realize this, let p_ (t) be the fraction of network links connecting agents 
with different opinions. The probability that any agent changes its opinion at 
a given step is n\{t) = and the probability that a connection is broken 

is n 2 (t) = £>„(t)(l — Pi)p 2 - The sum n(t) = ni(t) + n 2 {t) is the probability per 
evolution step that any change takes place, and thus fixes an overall evolution 
time scale. If p\ and p 2 vary in such a way that the ratios q = iii(t)/ii(t) and 
1 — q — 7i 2 (t)/n(t) are kept constant, such overall time scale will change, but 
the relative frequency of the two processes will be the same, giving rise to 
statistically equivalent final states. Thus, the only independent combination 
of the probabilities pi and p 2 relevant to the determination of the final state 
is 

q=?M = * . (1) 

7r(t) Pi + {1-Pi)p 2 V ; 



The behaviour at the two extreme values of q is immediately assessed. For 
q = (i.e. pi — 0), where no opinion changes take place, the final structure 
consists of two mutually disconnected communities with similar sizes and op- 
posite opinions. Internally, each of them stays fully connected, so that the total 
number of remaining connected pairs R is close to 2(N/2)(N/2 — 1)/2 « N 2 /4. 
For q = 1 (i.e. p\ = 1), on the other hand, interacting agents with initially 
opposite opinions always reach consensus, so that no interaction links are bro- 
ken. At the final state, all agents share the same opinion and the network is 
still fully connected, with its N(N - l)/2 « N 2 /2 links intact. 

These two limits suggest that the fraction of remaining connected pairs, r = 
2R/N(N — 1), provides a first quantitative characterization of the final state. 
Figure 1 shows r as a function of q, for systems of various sizes N. For each 
value of q, the fraction r has been obtained as the average over 500 to 5000 
realizations of the initial condition and the evolution. Rather unexpectedly, we 
find that r reaches a minimum for an intermediate value of q. The position g m j n 
and depth r min of this minimum depend on N, and seem to tend to zero as the 
population grows, as shown by the plot in insert (b). The network connectivity 
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Fig. 1. The fraction r of remaining links as a function of the parameter q, for N = 20 
(x), 50 (o), 100 (full line, and 500 (•). The dashed line represent the analytical 
approximation for large q discussed in the text. Insert (a): Close-up of the same 
data for small q. The dashed lines represent the small-g approximation discussed in 
the text. Insert (b): The position q mm (o) and the depth r m - m (•) of the minimum 
of r as functions of N. Straight lines are least-square fittings with slopes —0.28 (full 
line) and —0.75 (dashed line). 



at the minimum is considerably lower than at the extreme values q = and 
1. For N = 100, for instance, we have r min m 0.13, which implies an average 
connectivity of about 13 links per site. In this intermediate regime, thus, the 
network is poorly connected and the population structure correspondingly 
degraded. 

The minimum at g min defines two regimes in the parameter q. To the left (q < 
Qmm), as expected from the behaviour for q = 0, the population becomes split 
into two separate communities with similar sizes and opposite opinions. As q 
grows towards g min , the internal connectivity of these communities decreases. 
Simultaneously, several small separate clusters, each of them containing just 
a few agents (typically, less than 5 for A^ = 500), are also found at the final 
state. 



To the right of the minimum, we still have realizations where the population 
splits into disconnected communities. The two largest communities always 
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Fig. 2. The fraction / of realizations ending in a single community (•), and the 
contribution fro of those realizations to the fraction of remaining links (o), where ro 
is the average number of remaining links in a single-community final state. Curves 
ar Bezier splines, added for clarity. 

have opposite opinions, but their size is much more variable than for q < 
g m i n . Moreover, the number and size of small communities grow. At the same 
time, it becomes increasingly frequent to find realizations where the whole 
population stays connected into a single community. In these realizations, the 
final network is not fully connected, but its connectivity is significantly larger 
than in the cases where separate communities are formed. As a matter of fact, 
it is the contribution of these realizations which determines the growth of the 
fraction of connected pairs r for q > g min . As q keeps growing, the fraction 
of realizations with a single-community final state increases and, eventually, 
all realizations end on such states, with a homogeneous opinion all over the 
system. Figure 2 shows the fraction / of realizations where the population 
reaches a single-community final state as a function of q for N = 500, and the 
contribution of those realizations to the fraction of remaining links. 

In order to compare the regimes at both sides of the minimum, we have chosen 
to study in more detail the final structure of a population of size N = 500 
for two values of q where the fraction of remaining connected pairs r attains 
similar levels, namely, q = 5 x 10~ 3 and 0.3 (r 0.27; see Fig. 1). Figure 3 
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shows, for the two chosen values of q, the size of the second largest community 
(N 2 ) as a function of the size of the largest community (Ni), for 10 3 realiza- 
tions. Generally, N 2 ^ 0, except in realizations where the final state consists 
of a single community, where Ni = N and N 2 vanishes. The straight full line 
stands for the function N 2 — N — N±, so that the vertical distance from each 
dot to the line represents the size Nq — N — Ni — N 2 of the population not 
included in the two largest communities for the corresponding realization. The 
dashed line, in turn, corresponds to N 2 = (N — AT 1 )/2. Dots below this line 
represent realizations where the second largest community is smaller than N . 

For q < g m i n , we find that the size of the largest community reaches, at most, 
Ni ~ 340 ~ 0.7 N. Moreover, most dots lie over the straight line, which 
implies that essentially all agents are in one of the two largest communities 
(Nq pa 0). For q > g min , on the other hand, the largest community can reach 
the maximum size Ni = 500 = N, indicating that some realizations already 
correspond to single-community final states. Now, all sizes between Ni < N/2 
and N seem to be possible. Except in the close proximity of Ni = N, where 
the largest community comprises practically the whole population, dots are 
sensibly below the full line. The population N not belonging to the two largest 
communities is typically around 0.1N. Moreover, for most realizations where 
Ni > 400 = 0.87V, this population is above the size of the second largest 
community. In this situation, we can no longer properly speak of splitting 
into two main communities. The final structure is in fact closer to the single- 
community state, with the addition of several small separated communities. 

The distribution of internal connectivities in the resulting communities is also 
strongly dependent on the parameter q. Figure 4 shows the number of links 
Pi within each community as a function of the community size Ni, for the 
two largest communities in each realization (i = 1,2). Two different symbols 
identify the largest (o) and the second largest (•) communities. A third symbol 
(■) is used for realizations with a single-community final state (Ni = N). 
The full curve represents the function P^ 2 = Ni i2 (Ni >2 — l)/2, corresponding 
to fully connected networks. The dashed curve stands for Pi j2 = Ni y2 — 1, the 
minimum number of links in a connected network. 

For q < g m i n , the number of remaining links in each of the two largest com- 
munities exhibits a limited dispersion between realizations -just like their size 
N\ >2 . As a result, all the realizations are represented by dots in a rather com- 
pact cloud. Within each community, roughly one half of the total possible 
connections are actually present, P lj2 pa N lj2 (N 1:2 — l)/4, so that the aver- 
age connectivity per site is about N 1>2 /2. In this situation, thus, the internal 
connectivity of both communities is quite high. 

For q > q m i n , in contrast, the number of remaining links is broadly distributed. 
For practically all sizes, the number of links of second-largest communities 
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Fig. 3. The size N 2 of the second largest community as a function of the size Ni 
of the largest community, for two values of the control parameter q, in a series of 
10 3 realizations of a system with TV = 500 agents. Each dot corresponds to a single 
realization. The straight and the dashed lines represent, respectively, the functions 
N 2 = N - Ni and N 2 = (N - iVi)/2. 

is just above the minimum, P 12 = A^ 12 — 1, indicating a barely connected 
structure. The average connectivity per site within these communities is, in 
fact, close to one. The number of links in the largest communities, on the 
other hand, are always considerably above the minimum. However, even the 
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Fig. 4. The number of remaining links p within the largest (o) and the second 
largest (•) community (i = 1,2, respectively), as a function of the community size 
Ni, for two values of q, in a series of 10 3 realizations of a system of N = 500 agents. 
Realization with single-community final states are indicated by a different symbol 



(■). The full curve in each plot represents the function Pi^ 
dashed curve is Pi 2 = N± 2 — 1. 



^1,2(^1,2 -l)/2; the 



8 



better connected among these communities are still far from the situation 
of full connectivity. For realizations where Ni < N/2, P\ reaches at most 
10% of its maximum value. The number of links almost attains the maximum 
in some of the realizations where the population stays connected in a single 
community. These realizations, in fact, fill the gap between the best connected 
communities in the case of split populations and fully connected networks. 
Interestingly, while the above description emphasizes the differences between 
the internal structure of second- largest communities, largest communities, and 
single-community populations, the fact that all the dots in the lower panel of 
Fig. 4 seem to lie over a smooth curve suggest that those different structures 
belong to a sole class. The internal connectivity varies continuously as the 
community size changes. 

Obviously, the coevolution of the distribution of opinions and the interac- 
tion network creates strong correlations between the individual states of those 
agent pairs that remain connected. These correlations make an accurate ana- 
lytical treatment of our model particularly difficult. Some of the features found 
in the numerical results, however, admit an analytical explanation under suit- 
able hypotheses. Let us begin by the regime of large q where, as we have seen, 
most realizations end in a single-community state. The evolution towards a 
homogeneous opinion all over the population is fully driven by fluctuations. 
In fact, the number N + of agents with opinion +1, for instance, performs 
a random walk, N + — > iV + ± 1, each step an opinion changes. Starting at 
N + fa N/2, the process ends when N + reaches either or N. From standard 
results on first-passage time distributions in random walks [12], we find that 
the average number of active steps needed to reach one of these extreme values 
is S — N 2 /4px. An active step is here defined as an evolution step where any 
change, either the flip of an individual opinion or the deletion of a network 
link, occurs. In an active step, a link is deleted with probability (1 — p\)p-i- If 
we assume that this latter event is not correlated with opinion changes, the 
total number of links deleted in S active steps is N 2 (l — Pi)p2/4pi- Taking 
into account that the initial number of links is, approximately, N 2 /2, we have 



N 2 /2 



N 2 N 2 (l- Pl )p 2 
2 4 Pl 



Note that this estimation for the fraction of remaining links is independent of 
N. It is represented by the dashed curve in Fig. 1. We see that the agreement 
with numerical results is excellent for large q. While neglecting correlations 
between link deletion and opinion changes is not strictly correct -as a matter of 
fact, links are deleted precisely in those steps where opinions do not change- for 
q fa 1 deletion events are so rare that such correlations can hardly accumulate 
into a relevant effect before the final state is reached. 

For small q, on the other hand, the final state always corresponds to essentially 
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two large communities with opposite opinions. Since, in this limit, pi is small, 
most links connecting agents with different opinions are deleted before any 
significant number of opinion flips occurs. At any time, thus, essentially all 
links connect agents with the same opinion, except just after an opinion flip: 
when an agent's opinion changes, most of its links are now connections with 
agents with the opposite opinion. The total number P + of links connecting 
agents with the same opinion remains constant until an opinion flip takes 
place. In such event, the decrease in P + is given by the number of links of the 
agent whose opinion is changing, which can considerably vary among agents. 
As a kind of mean-field approach, we assume that this number is proportional 
to the average per-agent number of links connecting agents with the same 
opinion, 2P+/N, so that P + = —2ap\P + /N. The derivative is performed with 
respect to a time variable whose units are active steps, and the heuristic factor 
a represents the proportionality assumed above. The solution to this equation 
is 

P +(s) = ^exp(-?fl S ), (3) 



where s is the number of elapsed active steps. Moreover, as discussed above, 
the total number of links after s active steps is 

N 2 

P(s) = — -(1- Pl )p 2 s. (4) 



The evolution ends when P+(s) = P(s), so that the fraction r of remaining 
links is determined by the equation 



r = - exp 



—aN- 



Q 



1-q 



:i-r) 



(5) 



Insert (a) of Fig. 1 shows fittings of the numerical results for small q, obtained 
from this equation with a = 0.4. The fact that a does not depend on the 
size N implies that r depends on q and N through their product qN. This 
explains, at least qualitatively, that the minimum in r shifts to smaller q as iV 
grows. This shift is a direct consequence of the different dependence on N of 
the two mechanisms which drive the system for small q. While at each active 
step only one link can be deleted, the number of links connecting agents with 
equal or different opinions can change in a quantity of order N. 

In summary, we have shown that our simple model of coevolution for a pop- 
ulation of agents and the underlying interaction network gives rise to an in- 
teresting variety of population structures. The evolution rules represent the 
spreading of a bivaluate opinion on a network whose links can break when 
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agents do not succeed at reaching an opinion agreement. The final population 
structure consists, typically, of a set of separate communities, each of them 
containing agents with the same opinion. The resulting structures can be di- 
vided into three main classes: (i) two internally well connected communities 
with similar sizes and opposite opinions; (ii) a single community containing 
all the population; and (iii) a well connected community with, typically, more 
than half the population, accompanied by a set of poorly connected smaller 
communities. In this latter class, the two largest communities always have 
opposite opinions. 

While a clear boundary between the three classes cannot be unambiguously 
drawn, class (i) on the one hand, and (ii) and (iii) on the other, are observed 
within different ranges of the control parameter q 6 [0, 1]. The cross-over re- 
gion between these two ranges moves towards q = when the population size 
iV grows, so that the small-g regime would disappear for N — >■ oo. This ef- 
fect can be avoided if, when two agents fail to reach an opinion agreement, 
not only their mutual connection but a given fraction of their links to other 
agents is also broken. In this way, deletion of links and opinion change would 
have the same dependence on the system size. Preliminary numerical results 
of this extension of the model have already been obtained, and will be pub- 
lished in a forthcoming paper. Another extension, in the direction of making 
the model more realistic, would allow for the creation of links between discon- 
nected agents. This process would drive the population to a dynamic asymp- 
totic state, independent of the initial structure, where communities could form, 
aggregate, exchange agents, and disappear, as known to happen in real social 
systems. 
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